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Andreev reflection at the interface between a half-metallic ferromagnet and a spin-singlet su- 
perconductor is possible only if it is accompanied by a spin flip. Here we calculate the Andreev 
reflection amplitudes for the case that the spin flip originates from a spatially non-uniform magneti- 
zation direction in the half metal. We calculate both the microscopic Andreev reflection amplitude 
for a single reflection event and an effective Andreev reflection amplitude describing the effect of 
multiple Andreev reflections in a ballistic thin film geometry. It is shown that the angle and energy 
dependence of the Andreev reflection amplitude strongly depends on the orientation of the gradient 
of the magnetization with respect to the interface. Establishing a connection between the scattering 
approach employed here and earlier work that employs the quasiclassical formalism, we connect 
the symmetry properties of the Andreev reflection amplitudes to the symmetry properties of the 
anomalous Green function in the half metal. 

PACS numbers: 74.45.+c,74.78.Na,75.70.Cn 



I. INTRODUCTION 

Superconductors extend their order into adjacent 
normal- metals via the mechanism of Andreev reflection: 1 
An electron incident on the superconductor-normal- 
metal interface is phase-coherently reflected as a hole, 
and vice versa. 2 Half-metallic ferromagnets (half metals) 
are normal metals that support quasiparticle excitations 
of only one spin orientation. If superconducting corre- 
lations are to extend from a spin-singlet superconduc- 
tor into a half-metallic ferromagnet, the azimuthal spin- 
rotation symmetry around the magnetization axis needs 
to be broken near the superconductor-half-metal inter- 
face, so that electrons can be Andreev reflected into holes 
in the same spin band. 

In standard ferromagnets (which have majority as well 
as minority carriers), such a combination of spin-flip 
and Andreev reflection was shown theoretically 3-6 and 
experimentally 7-10 to lead to a superconductor proximity 
effect with a range comparable to that in normal metals. 
In the literature, magnetic domain walls, 3-5 ' 11-14 helical 
structures intrinsic to the ferromagnetic material, 1 
artificially structured multilayers with noncollinear mag- 
netization directions in different layers, 18-22 and a pro- 
cessing magnetization direction 23 were addressed as pos- 
sible microscopic origins of the broken spin-rotation sym- 
metry. Although the possibility of spin-flip Andreev re- 
flection (with similar origins) also exists for half-metal- 
superconductor interfaces, 20 ' 24-29 and a sizeable Joseph- 
son current has been observed in junctions involving 
a half-metallic ferromagnet, 30,31 the Andreev reflection 
mechanism at half-metal-superconductor interfaces was 
found to be more delicate than in the case of standard 
ferromagnets. 29 

The underlying reason for the differences between a 
ferromagnet-superconductor (FS) interface and a half- 



metal-superconductor (HS) interface is that, if no orbital 
symmetries are broken, the latter admits a description in 
terms of an effective 2x2 scattering matrix, whereas 
the former requires at least a 4 x 4 scattering matrix, to 
account for the spin degree of freedom. Mathematical 
constraints on 2 x 2 matrices following from current con- 
servation and particle-hole symmetry then force the An- 
dreev reflection amplitude rh e (e) of an HS junction to be 
generically zero at the Fermi level e = 0, whereas there is 
no such strong restriction on the Andreev reflection am- 
plitudes for an FS junction. As shown in Ref. 29, the fact 
that rh o (0) = for an HS junction has little consequences 
for the strength of the proximity effect at distances be- 
low the superconducting coherence length £sj but leads 
to a suppression of the proximity effect at larger dis- 
tances in comparison to a normal-metal-superconductor 
(NS) junction. 29 For example, the zero-temperature crit- 
ical current in a ballistic SHS junction of length L 3> £s 
is a factor ~ (L/£s) 2 smaller than in an otherwise com- 
parable SFS junction. 29 (The experiments of Refs. 30,31, 
which observe a sizeable Josephson effect in SHS junc- 
tions, have L ~ £s, so that they need not be affected by 
this effect.) 

If orbital symmetries are broken, the mathematical 
constraints leading to the condition r^cfe) = at e = 
are no longer operative. Examples of broken orbital sym- 
metries that allow for a finite Andreev reflection ampli- 
tude rho (e) at e = are a magnetization gradient parallel 
to the interface, spin-orbit coupling in the superconduc- 
tor — which both break inversion symmetry with respect 
to the interface normal — , or impurity scattering. An- 
dreev reflection in the presence of a magnetization gra- 
dient parallel to the interface, which is relevant for a do- 
main wall in the half metal, was previously considered by 
us in Ref. 32 in the limiting case of a half metal with infi- 
nite wavefunction decay rate kj, of the minority carriers. 
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This situation will be analyzed in more detail here. We 
lift simplifying assumptions of Ref. 32, and compare the 
generalized result to the case when the magnetization 
varies perpendicular to the interface, where no orbital 
symmetry is broken. The other two scenarios for break- 
ing orbital symmetries will be analyzed elsewhere. 33,34 

The scenario of a magnetization gradient parallel to 
the HS interface is particularly relevant for a lateral 
superconductor-half- metal contact, for which the super- 
conducting contact is deposited on top of a half-metallic 
film. The experiments of Refs. 30,31 were performed in 
such a lateral geometry. A lateral contact, with a do- 
main wall in the half metal under the contact area, is 
shown schematically in Fig. 1. For a phase-coherent lat- 
eral contact, it is convenient to define an "effective" An- 
dreev reflection amplitude r£f that represents the effect 
of multiple Andreev reflections at the HS interface for 
a quasiparticle moving in the thin half-metallic film, in- 
cident on a domain wall (or any other region with non- 
uniform magnetization) from a region with uniform mag- 
netization, see Fig. lc. (The standard definition of the 
Andreev reflection amplitude r^ c is for single incidence of 
a quasiparticle onto the superconductor interface. This 
is the situation shown in Fig. lb.) Because one is in- 
terested in currents flowing parallel to the interface in a 
lateral contact (see Fig. lc), it is the effective amplitude 
r^f, not The, which occurs in a calculation of, e.g., the 
subgap conductance or the Josephson current. 32 

In this article, we calculate both the microscopic An- 
dreev reflection amplitudes and the effective ampli- 
tude r£g . We consider an HS interface with a low normal- 
state transmission probability, caused by the presence of 
a tunnel barrier or by the mismatch of Fermi velocities on 
both sides of the junction. We do not consider the effect 
of impurity scattering. For the calculation of r^ c this is 
not a serious shortcoming as Andreev reflection is local, 
and only impurities in the immediate vicinity of the in- 
terface would play a role (up to distances of the order of 
a Fermi wavelength) . For the effective Andreev reflection 
amplitude r^f in the thin-film geometry, the omission of 
impurity scattering limits the applicability of our results 
to situations in which the mean free path is larger than 
the domain wall size Id- In addition to the scenario of a 
magnetization gradient parallel to the HS interface, we 
also consider the case that the magnetization gradient is 
perpendicular to the HS surface as is relevant, e.g., if the 
half metal has different boundary and bulk anisotropies. 

The calculation of the microscopic, single-reflection 
amplitude rh c will be presented first. Hereto, we first 
review the restrictions of symmetry on Andreev reflec- 
tion at a HS interface in Sec. II. We then introduce the 
Hamiltonian and the wavefunctions for the case of a uni- 
form magnetization in Sec. III. The calculation of the 
Andreev reflection amplitudes in the presence of a non- 
uniform magnetization then takes place in Sec. IV. The 
calculation of the effective Andreev reflection amplitude 
r ho f° r a superconductor placed on top of a thin half- 
metallic film is then presented in Sec. V. Finally, in Sec. 
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FIG. 1: (Color online) (a) Lateral contact between a thin 
half-metallic film (H) and a superconductor (S). Andreev re- 
flection at the HS interface is possible at those interface po- 
sitions where the magnetization direction (indicated by thin 
black arrows) is not uniform, such as a domain wall. The 
situation shown in the figure is generic, since domain walls 
at the HS interface are bound to occur if the HS contact is 
larger than the domain size in the half-metallic film. In this 
geometry, there are two possible ways to define the Andreev 
reflection amplitude: The amplitude rhc for a single Andreev 
reflection (thick red arrows), appropriate for the interface be- 
tween a semi- infinite half metal and a superconductor (b), and 
the effective amplitude rjf representing the combined effect 
of multiple Andreev reflections in the same region of non- 
uniform magnetization for quasiparticles moving in the half 
metallic film (c). In parts (b) and (c), the black arrows in 
the top panel show the spatial variation of the magnetization 
direction, whereas red arrows in the bottom panel show an 
example of a trajectory for an incoming electron (solid) and 
an Andreev reflected hole (dashed). 



VI, we discuss the relation of the scattering approach 
used in this article and the Green function approach used 
in most of the literature. In particular, we consider the 
anomalous Green function and its frequency dependence. 
We conclude in Sec. VII. The appendix contains certain 
details of the calculations not presented in the main text. 



II. CONSTRAINTS IMPOSED BY UNITARITY 
AND PARTICLE-HOLE DEGENERACY 

Although the considerations of this section are com- 
pletely general, for definiteness we choose coordinates 
such that the half-metal-superconductor interface is the 
plane z = 0. In this section and in the following two 
sections, the half metal occupies the half space z < 0, 
and the superconductor, which is taken to be of s-wave, 
spin-singlet type, occupies the half space z > 0. This 
setup is shown schematically in Fig. 2. 
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FIG. 2: (Color online) Illustration of the scattering 
setup. Coordinates are chosen, such that the half-metal- 
superconductor interface is the plane 2 = 0. Electron-like 
quasiparticles incident on the HS interface are either reflected 
as an electron (normal reflection), or as a hole (Andreev re- 
flection). At a translationally invariant interface, the projec- 
tion k|| of the wavevector on the interface is conserved upon 
reflection. 



Quasiparticles in the half metal are labeled by their 
wavevector k|| = k x e x + k y e y parallel to the HS interface 
(e^ and e y are unit vectors in the x and y directions, re- 
spectively) and by their excitation energy e. We assume 
periodic boundary conditions in the x and y directions, 
so that the wavevectors ky are discrete. 

At excitation energies e below the superconduct- 
ing gap A, quasiparticles incident on the half- metal- 
superconductor interface from the half-metallic side will 
be reflected back into the half metal. This reflection 
can be either normal reflection or Andreev reflection, for 
which electron-like quasiparticles are reflected as holes 
and vice versa. The reflection process is described by a 
scattering matrix <S(kJ,,k|i;e), which takes the form 



, , . ee (k||,ky;e) r eh (kj p ky; e) 



(1) 



where the subscripts e and h refer to electron-like and 
hole-like states. The scattering amplitudes r c h and The 
describe Andreev reflection processes. 

The scattering matrix <5(kji , ky ; e) satisfies two con- 
straints: Unitarity and particle-hole symmetry. The lat- 
ter condition reads 



1 




(2) 



Sm,k f e)=(° 1 J)s(-kf|,-k|,;- e r(2 

The combination of unitarity and particle-hole degener- 
acy severely restricts the form of the scattering matrix if 
there is translation invariance along the interface, which 
implies 

S( k ||> k ||;e) = 5(ky;e)5 k j |ik|| . (3) 
Combination of Eqs. (2) and (3) gives 

r ee (k||,£) = r hh (-k||, -e)*, 

r e h(k||,e) = r hc (-k||,-e)*. (4) 

If the scattering problem is also invariant for a n rotation 
around the interface normal, so that 



one finds, upon combining all symmetry properties, 

r cc (k||;0)r ch (k||;0) =0 (6) 

for the scattering matrix at the Fermi level e = 0. Since 
r cc =/= 0, except for a special choice of parameters, this 
implies that generically one must have 29 



r eh (k||;0) =0. 



(7) 



These very general considerations pose a severe re- 
striction on the magnitude and the spatial extension 
of the proximity effect in half metals that is absent 
in ferromagnet-superconductor junctions with otherwise 
comparable characteristics. A nonzero Andreev reflec- 
tion amplitude for a half-metal-superconductor junction 
can be obtained only by fine-tuning device parameters 
such that the normal reflection amplitude becomes zero, 
or by invoking processes that break the symmetries lead- 
ing to Eq. (7). The former scenario was discussed in 
Ref. 29 and will not be addressed here. Examples of 
symmetry-breaking processes that result in a nonzero 
Andreev reflection amplitude are: lifting of particle-hole 
degeneracy by a finite excitation energy e, 29,35 breaking 
of the rotation symmetry around the interface normal, 32 
breaking of the translation symmetry along the interface, 
or the breaking of phase coherence. 36 A domain wall for 
which the magnetization direction varies in a direction 
parallel to the interface is an example of a perturbation 
that breaks the rotation symmetry. 32 However, a thin in- 
terface layer of different magnetic orientation than the in- 
terior of the half metal (which is a model of a "spin-active 
interface" 24 ~ 27,37 ) does not lift the constraints leading to 
Eq. (7). 29 ' 35 The role of variations in the magnetization 
direction will be considered in more detail in Sec. IV be- 
low. 

A finite excitation energy e lifts the particle-hole de- 
generacy, and the Andreev reflection amplitude r e h be- 
comes nonzero. If no other symmetries are broken, the 
order of magnitude of the Andreev reflection amplitudes 
at finite e can be estimated as 



kch(e)| 



min(A/r, E^ 



Teh.FSl 



(8) 



5(kj| ) k|,;e)=5(-kf„-k,|;e) ) 



(5) 



where E% is the Thouless energy of the interface layer 
where the singlet-triplet conversion takes place, r is the 
transparency of the superconductor interface, A the su- 
perconducting gap, and r e h, fs the Andreev reflection am- 
plitude of a ferromagnet-superconductor amplitude of 
otherwise comparable characteristics. The first energy 
scale in the denominator comes about because electrons 
and holes scattering off a normal-metal-superconductor 
interface of transparency r at finite excitation energy 
e experience an additional phase difference ~ ±er/A, 
which lifts the electron-hole degeneracy. 38 The second en- 
ergy scale in the denominator appears from phase differ- 
ences acquired in the interface layer. The typical thick- 
ness of this interface layer is of the order of minority 
decay length £, which implies that E^ is of the order of 
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the Fermi energy. For tunneling interfaces one always 
has |r e h(£)| <C l^oh.Fsl and we conclude that the break- 
ing of electron hole symmetry by finite excitation energies 
is not an efficient route towards sizeable Andreev reflec- 
tion in that case. The suppression of Andreev reflection 
in half-metal-superconductor junctions (as compared to 
ferromagnet-superconductor junctions) is absent only for 
transparent interfaces and excitation energies of order A. 

The ^-dependence of r c h not only determines the 
conductance through the half-metal-superconductor in- 
terface at finite bias, 29 ' 35 it also sets the scale for 
the Josephson effect in a superconductor-half-metal- 
superconductor junction. 39,40 If the Thouless energy El 
of a Josephson junction of length L is large in compari- 
son to A ("short junction limit"), the Josephson current 
/ is carried by quasiparticle states with energies up to 
A. In this limit, the symmetry considerations that sup- 
pressed Andreev reflection at e — do not affect the or- 
der of magnitude of /, and one concludes that otherwise 
comparable superconductor-half- metal-superconductor 
and superconductor-ferromagnet-superconductor junc- 
tions have comparable Josephson currents. 29 If, how- 
ever, El <C A ("long junction limit"), only quasipar- 
ticle states with energy below El contribute to /, so 
that I is significantly suppressed compared to the Joseph- 
son current in comparable superconductor-ferromagnet- 
superconductor junctions. 

In the remainder of this article, we present explicit 
model calculations of the Andreev reflection amplitudes 
for the case that singlet-triplet conversion is mediated by 
a spatially non-uniform magnetization in the half metal, 
as it appears, e.g., in a domain wall. 



III. HAMILTONIAN AND SCATTERING 
STATES 

A. Bogoliubov-de Gennes Hamiltonian 

Quasiparticle excitations near the HS interface are de- 
scribed by the Bogoliubov-de Gennes equation 41 

where the four-component spinor 

*(r) = ( Ut (r),^(r),« T (r), n (r)) T (10) 

consists of wavefunctions u a (r) for the electron and 
v a (r) for the hole degrees of freedom. The supercon- 
ducting order parameter A(r)e l< ^ is nonzero only in the 
superconductor. We will take A(r) = A9(z), where 
®{z) = 1 if z > and otherwise. This step function 
model is appropriate for tunneling interfaces of s-wave 
superconductors. 42 

For the single-particle Hamiltonian, we take the sim- 
plest model that contains the essential features of the 



half-metal-superconductor interface, 

H = -ft 2 V^-yV-^^(z)P ff (r) + ^(z), 

(11) 

where 

m(z) = ( mH l Z < °> (12) 
w mg if z > 0, v ' 

with ran and ms being the effective masses for the half 
metal and the superconductor, respectively, 

^-{ns ifz>0, (13) 

with a =t,i and the potentials fj,^, yUnj., and repre- 
senting the combined effect of the chemical potential and 
band offsets for the majority and minority electrons in 
the half metal and for the superconductor, respectively, 
and where w sets the strength of a delta-function poten- 
tial barrier at the interface. The operators 

A = ^ + ^m(r) • <x, n = i-im(r).<x (14) 

project onto the majority and minority components, re- 
spectively, where m(r) is a unit vector pointing along the 
magnetization direction in the half metal. 

The potentials £iHt> Mhj., and us are chosen such that 
A'Hti Ms > 0, and /lihj, < 0. As a result, majority states in 
the half metal and in the normal state of the supercon- 
ductor are propagating states, with Fermi wavenumbcrs 



fef = ^V 2m HA*Hti k s = -^msMS, (15) 

respectively. The corresponding Fermi velocities are Vf = 
hk^/mn and vs = hks/ms, respectively. Minority states 
in the half metal are evanescent with wavefunction decay 
rate 



ki = -y/2m H \mi\- (I 6 ) 

(The wavefunction decay rate is the inverse of the 
wavefunction decay length £ used in the previous sec- 
tion.) The strength w of the (5-function potential is cho- 
sen such that the transmission probability of the inter- 
face is much smaller than unity. It is in this limit only, 
that the step-function model for the superconducting or- 
der parameter A used in Eq. (9) is valid. 42 We use the 
Andreev approximation A <C throughout our calcu- 
lation. 



B. Scattering states for A = 

In order to introduce the relevant notation, we first 
consider solutions of the Bogoliubov-de Gennes equation 
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(9) in the normal state (i.e., with A = 0), for a spatially 
uniform magnetization direction m = e 3 , and at e = 
0. In this case, solutions of the Bogoliubov-de Gcnnes 
equation (9) can be written as a product 

*(r) = e* k ir r v]/ k|| ( z ), (17) 

where ky = k x e x + k y e y . For z < 0, the spinor wavefunc- 
tion 'Ski, has the general form 



*k„(*) - 



Chfe" 



ik ^ z + c^e ik ^ z 



The amplitudes r, r' , and t are majority electron re- 
flection and transmission amplitudes of the half-metal- 
superconductor interface (with the superconductor in the 
normal state); the amplitude r'^ is the minority elec- 
tron reflection amplitude. The coefficient t± parameter- 
izes the evanescent wave amplitude for minority electrons 
in the half metal. (There is no transmission amplitude 
for minority electrons.) In terms of these amplitudes, 
the assumption of a tunneling interface translates to \t\, 

\h\ « i. 



/Viz 






Chie K ^ z 



where 



Hz 



and 



(18) 



(19) 



(20) 



In the superconductor, for z > 0, the general form of the 
spinor wavefunction is 



/ d cf e- lks * z + d' ct e iks * z \ 



*k„(*) 



'VSz 



ik$ z z 



d ei e~ u 
d h1 -e 

V d H e ik *.* + d' H e- ik ^ z j 



7 „— ikszZ 



d'ht e 



where 



ksz 



VSz 



Hksz/m s . 



(21) 



(22) 



The amplitudes appearing in the above equations are re- 
lated as 



t r' 
r* t* 
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d' ei 




Chi 










(23) 



with 



2iw + Vf z + v Sz ' 

-1 + t\Jvf z /vSz, 



-l + t^Jv Sz /v tz , 

-1 + t l^J VSz /Viz ■ 
Zy/VlzVSz 

2iw + iv iz + v Sz ' 



(24) 



C. Scattering states for uniform magnetization 



We now use the notation established in the previous 
subsection to construct retarded and advanced scatter- 
ing states for the half-metal-superconductor interface at 
finite excitation energy e. The scattering states will be 
used for the perturbation-theory calculation of the An- 
dreev reflection amplitudes for a non-uniform magneti- 
zation in the next section. 

As before, we consider a spatially uniform magnetiza- 
tion direction m = e3. For each wavevector k||, there 
is an electron-like and a hole-like scattering state, which 
we label |k|| , e) R and |k|| , h) R for the retarded states and 
|k||,e) A and |k||,h) A for the advanced states. A retarded 
scattering state is called "electron-like" or "hole-like" if 
the incoming part is electron-like or hole-like, respec- 
tively, whereas one considers the outgoing part of the 
wavefunction for an advanced scattering state. In gen- 
eral, the outgoing part of a retarded scattering state will 
be of mixed electron/hole type, as well as the incom- 
ing part of an advanced scattering state. For the states 
constructed below, however, there is no mixing between 
electron-like and hole-like parts in the propagating com- 
ponents, since there is no Andreev reflection at an HS in- 
terface if the magnetization is uniform. Only the evanes- 
cent parts of the scattering states will be of mixed elec- 
tron/hole type. 

The scattering states are constructed from the general 
form of the spinor wavefunction (z) at finite excita- 
tion energy e and with nonzero superconducting order 
parameter A, which reads 



/ c cf e lk ^^ z + c' ct e- tk ^^ z \ 



*k„ (z) 











/Vjz 





Cci ^ u (e)z 



e Klz{-£)z 



(25) 
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for z < 0, and 



*k„(*) - 



= ifcs2Z-K Sz (e)z 
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f p -"?(e)-*^ 
4- 



+ ■ 



-iks z z~Ks z {e)z 

V 7 ^ 



( 



(If 



-die +ir >&- i ' 1 ' 



for z > 0. Here we denned 



s/hv tz , 



fc Tz (e) = hf z 

77(e) = arccos(e/A), 
Ks ,(e) = (VA 2 ^ 2 )/fo S! 



(26) 



(27) 
(28) 
(29) 

(30) 



(In Eqs. (27), (28), and (30) we kept the leading terms in 
a small-e expansion only.) Solution of the Bogoliubov-de 
Gennes equation (9) in the Andreev approximation then 
gives the following relations between the coefficients 
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d' 



r + 



e 2ir)(e) r ^ _ r i 



Cef, 



tt\.e 



ir}{e)—i<f> 



e 2i V (s) r > _ r / Ce t' 
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(31) 



With the help of the above wavefunctions, we construct 
the retarded scattering states |ky , e) R and |k|| , h) R as the 
state with wavefunction 



(r|k| h e) R 
(r|k|,,h) R 



1 



*k„,e(«)e' 



*k,,,h(z)e' 



ikii -r 



ki r r 



(32) 



where W x and W y are the dimensions of the HS inter- 
face in the x and j/ directions, the spinor wavefunction 
^k«,c(z) is given by Eqs. (25) and (26) with c e t = 1, 



Chi = 0, all other coefficients being determined by Eqs. 
(31), and the spinor wavefunction \&k||,h(z) is given by 
Eqs. (25) and (26) with c e t = 0, Chi = 1. Similarly, 
the advanced scattering states |k||,e) A and |kii,h) A are 
then defined as the states for which c' c ^ — 1, c h ^ = and 
Cg t = 0, c[, = 1, respectively, again with all other coeffi- 



cients determined by Eqs. (31). The factors 1/ yjW x W y 
in Eq. (32) ensure that retarded and advanced scattering 
states are normalized to unit incoming or outgoing flux 
in the z direction, respectively. The lengths W x and W y 
do not appear in the final expressions for the Andreev 
reflection amplitudes. These scattering states are at the 
basis of the calculation of the Andreev reflection ampli- 
tudes in the presence of a non-uniform magnetization, 
which is described in the next section. 



IV. ANDREEV REFLECTION IN THE 
PRESENCE OF A NON-UNIFORM 
MAGNETIZATION DIRECTION 

A. Slow variations of the magnetization direction 

A spatial variation of the magnetization direction in 
the half metal breaks the remaining symmetries in spin 
space and allows for Andreev reflection at the half-metal- 
superconductor interface. Here we consider a continu- 
ous and slow variation of the the magnetization direction 
m(r). An example of such a continuous change is a do- 
main wall, for which the net change of the magnetization 
angle is 7r, smeared out over a length much larger than 
the microscopic lengths kg , fcjT 1 , and k^ 1 . However, 
smaller rotation angles are possible, too, e.g. induced by 
strain at the interface due to lattice mismatches. 43 

To be specific, we choose a right-handed set of unit 
vectors ei, e 2 , and e 3 in spin space (which need not co- 
incide with the coordinates used for the orbital degrees of 
freedom) and consider a variation of the magnetization 
direction m of the form 

m(r) = (ei cos <j> m + e 2 sin 4> m ) sin 6 m (r) + e 3 cos m (r) . 

(33) 

Such variations of the magnetization direction are suf- 
ficient to model domain walls, but they do not allow 
for certain continuous changes of the magnetization at 
a fixed polar angle, as it occurs in helical magnets. (The 
full expressions for arbitrary variations of m are given in 
App. A.) We then employ a gauge transformation that 
rotates m to the e 3 -direction, 

H -+ W(r)t«W(r), U(v) = ( U ^ ^ ) , (34) 

with 

U(r) = e i0 m(m(r)xe 3 )-<T/2sin0 m ^ 

This gauge transformation adds a spin-dependent gauge 
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potential 



A(r) = 



ihuWu 

- (a 2 cos q 



)V0„ 



(36) 



to the Hamiltonian H 44 but it does not affect the sin- 
glet superconducting order parameter, since U' T ia2^U = 
ia^A. To lowest order in the rate of change of the angle 
9 m we then find that the perturbation V to the Hamilto- 
nian H reads 



(a) 



H 



71 7 7 7 7 7 

^ ^ 

yjS> /|V ^\ ^ 



H 



(b) 



1 4- 



FIG. 3: (Color online) Illustration of a magnetization gradi- 
ent perpendicular to the HS interface (a) and parallel to the 
HS interface (b). 



V — i(er 2 COS< 



<j\ sm ( 



h 2 H 2 
2m 2m 



(37) 



Since we take the length scale Zd for variations of the 
gradient of the magnetization angle 9 m to be large in 
comparison to the microscopic length scales kl , kg 1 
and k7 , we may neglect spatial variations of the per- 
turbation V in the direction parallel to the interface. In 
this approximation translation symmetry along the in- 
terface is preserved (after the gauge transformation) and 
the scattering matrix <5(k|| , k|| ; e) is diagonal, see Eq. (3). 
To lowest order in the rate of change of 6 m , the Andreev 
reflection amplitudes may then be calculated in pertur- 
bation theory. Using the scattering states defined in the 
previous section, one finds 



r hc (k||,e) = -^ A (k||,h, £ |V|k||,e, £ ) R , 
r e h(k||,e) 



h 

4 A (k||,e,e|V|k||,h, £ ) 



where 



V = 



V 

-V* 



(38) 



(39) 



We now present calculations of the Andreev reflection 
amplitudes for two special cases: Variation of the angle 
6 m in a direction perpendicular to the superconductor 
interface, as illustrated in Fig. 3a, and variation of 9 m 
in a direction parallel to the superconductor interface, 
which is illustrated in Fig. 3b. These two cases differ 
with respect to the symmetries discussed in Sec. II: The 
rotation symmetry around an axis normal to the interface 
is preserved in the former case, whereas it is broken in the 
latter case. We will see that this difference has profound 
consequences for the Andreev reflection amplitude. 

A variation perpendicular to the interface is rele- 
vant e.g., if the interface anisotropy at the half- metal- 
superconductor interface differs from the anisotropy in 
the bulk of the half metal. Calculating the Andreev re- 
flection amplitude, we find, to leading orders in \t\ and 



»"he(k||;£) 



■ £e - 4 (0-0 m ) A 



(V0 m • e,) 



8VA 2 



hk tz {v 2 z + v 2 z )^v^vTz 



(40) 



The dependence of this result on the interface parame- 
ters agrees with what was derived in Sec. II using gen- 
eral considerations. Note that in the limit — > oo, in 
which the minority carriers are completely expelled from 
the half metal, the Andreev reflection amplitude r^c van- 
ishes. The appearance of the azimuthal magnetization 
angle m as a shift of the superconducting phase <f> was 
found previously in Refs. 20,29 for the interface between 
a magnetic multilayer and a superconductor. 

A variation of the magnetization direction in which 
V0 m is parallel to the superconductor interface breaks 
the rotation symmetry around an axis normal to the in- 
terface, thus allowing, in principle, for a nonzero Andreev 
reflection amplitude at the Fermi level e = 0. 32 Here we 
elaborate on our previous calculation of this effect and 
generalize the results of Ref. 32 to the case of a finite 
minority wavefunction decay rate kj, in the half metal. 
Calculating the Andreev reflection amplitude according 
to Eq. (38), we then find, again to leading orders in \t\ 



and t 



4-1) 



nie(k||; £ ) 



_ (V0 m -k||) 



.(41) 



The first term in Eq. (41) comes from the overlap integral 
in Eq. (38) inside the superconductor, whereas the second 
term comes from the overlap integral in the half metal. 
The existence of a finite contribution to ri lc from inside 
the superconductor is responsible for the fact that 
remains nonzero in the limit kj, — > oo if 9 m varies parallel 
to the interface. In both cases, the amplitude for the 
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conversion of holes into electrons is given by 

r e h(k||;e) = r hc (-k|| ; -e)*, 

as discussed in Sec. II. For a general direction of V0 m , 
the Andreev reflection amplitudes are the sums of the 
contributions of Eqs. (40) and (41) above. 

We note that the presence of Andreev reflection at the 
Fermi energy for a domain wall with V(9 m parallel to 
the interface is accompanied by a nontrivial angle depen- 
dence of the Andreev reflection amplitudes r h and r^: 
If V0 m is parallel to the interface, rh c and r e h are even 
functions of e, but odd functions of k|| . On the other 
hand, if V# m is normal to the interface, r^c and r c h are 
odd functions of e, but even functions of k||. Similar 
behavior has been noticed previously for the anomalous 
Green functions. 27 ' 45 



B. Spin-active interfaces 

As a second example of a spatially varying magneti- 
zation direction, we now investigate a simplified model 
of a thin ferromagnetic or half-metallic layer located at 
the interface, the magnetization of which is misaligned 
with respect to the bulk of the half-metal. In this model 
of what is more generally referred to as a "spin-active 
interface" 24 - 27 > 37 W e take the magnetization direction to 
be the unit vector e3 in the entire half metal and consider 
a perturbation to the Hamiltonian H of the form 



V = hrh ■ &5(z), 



(42) 



where 



m 



(ei cos i 



e 2 sin </> m ) sin m + e 3 cos 9 m . (43) 



Spin-flip Andreev reflection at such "spin-active" HS in- 
terfaces has been considered previously in Refs. 27-29. 
Using Eq. (38) to calculate the Andreev reflection ampli- 
tude to first order in h and taking the limit of a tunneling 
interface, |ft^| 2 <C 1, we then find 



nic(k 



sm( 



2% z (A2- £ 2) 



(44) 



The proportionality to the square of the tunneling prob- 
ability is in agreement with the general considerations 
of Sec. II. The opposite limit of an ideal interface (\t\ — 
\t±\ = 1) was considered in Ref. 29. 



V. THIN FILM GEOMETRY 

A modification of the scattering problem arises in a 
lateral contact geometry as in Fig. la. For this geometry 
a formulation with scattering states describing quasipar- 
ticles moving in the plane of the half metallic film is more 
relevant than the formulation of the previous sections in 



> 



H t ^ 1, 



FIG. 4: (Color online) Illustration of the thin-film geometry. 
The coordinates are chosen such that the magnetization di- 
rection does not change in the y direction (perpendicular to 
the plane of the drawing). 



terms of scattering states of quasiparticles incident on 
the HS interface. 

The geometry we consider is shown in Fig. 4. As be- 
fore, the half-metal-superconductor interface is the plane 
z = 0. The superconductor occupies the half space 
z > 0, whereas the half metallic film is in the region 
—d < z < 0. The Hamiltonian is given by Eq. (9), and 
hard- wall boundary conditions are applied at z = — d. 
We choose our coordinates such, that (locally) the mag- 
netization direction does not depend on the coordinate y, 
so that the wavevector component k y is conserved, and 
consider the scenario that the magnetization direction m 
has a spatial non-uniformity only in a bounded region 
near x = 0. As in the previous sections, we consider the 
limit of a tunneling interface, and present our results to 
leading order in the transmission probability. 

We first construct the scattering states in the pres- 
ence of a uniform magnetization direction m = e,3 . The 
scattering states are normalized to unit flux in the x di- 
rection. There are electron-like and hole-like scattering 
states \k y n z s;e) and \k y n z s;h}, each labeled by the con- 
served wavevector component k y and by integers n z and 
s, where n z — 1,2,... represents the quantized transverse 
modes in the z direction and s = ±1 indicates scattering 
states propagating in the positive (s = 1) and negative 
(s = — 1) x direction. In the limit K±d 3> 1 the hard- wall 
boundary conditions at z = —d are inconsequential for 
the minority carriers, and the corresponding wavefunc- 
tions read, 



*fc y n^c(r) = -'< 



2v tx W y d 



*k,,,e(z)e 



ik x (s)sx-\-ik v y 



2v^ x W y d 



*k,,,h(2 ; )e 



—ik x (—e)sx-\-ik y y 



(45) 



where $]j ||jC (z) and ^k Ut h(z) are given in Eq. (32), 



h 2 



s 



e/hv tx , (46) 

and the allowed wavenumbers k^ z are determined from 
the condition that the wavefunctions vanish at z = — d, 
which implies the condition (compare with Eq. (31)) 



2k^ Z}B d + 7r — arg 



t 2 



(47) 



9 



for electron-like scattering states, and analogously for 
hole-like states 



2fc tZih d + 7T - arg r* + 



*2 



2nn z . (48) 



In the limit of a tunneling interface, one has r — > — 1 and 
\t\ <C 1, so that these conditions simplify to 

7TTI 

k u = -^+0{\t\), n z = l,2,... (49) 
and the wavefunctions of the scattering states become 



(r) 



v tx W y d 



e ik x (e)sx+ik yV gin "^ fz 



^k y n z s,h( r ) 



v tx W y d 



ik x {-e)sx+ik v y gin "z^ 2 



(50) 



For the purposes of our calculations higher-order terms 
in \t\ need to be kept in order to find a finite overlap 
between the different scattering states. In Eq. (45) we 
imposed periodic boundary conditions with period W y in 
the y direction in order to ensure a proper normalization 
of the scattering states. The normalization factors in Eq. 
(45) are valid in the limit of a tunneling interface only. 

We now consider the effect of a region in which the 
magnetization direction is not spatially uniform. In this 
case, scattering between the electron-like and hole-like 
quasiparticle states is possible. To lowest order in the 
rate of change of the magnetization direction, the An- 
dreev reflection amplitudes and for a quasiparti- 
cle incident on a region of nonuniform magnetization can 
be calculated in perturbation theory as 

rle (k y ;e;n' z s';n z s) = -^{k y n' z s';h\V\k y n z s;e), 

rl^(k y - : e- : n' z s';n z s) = --(k y n' z s';e\V\k y n z s;h), (51) 

where V is given in Eqs. (37) and (39) above. With these 
equations, the problem of calculating Andreev reflection 
coefficients is brought into a form similar to that of the 
previous section. 

In principle, the perturbation V depends on the coordi- 
nate z, even if the magnetization angle 9 m does not. [This 
can be seen, e.g., from Eq. (37).] Because of this, the 
scattering matrix is not diagonal in the transverse mode 
index n z . However, if the length scale I a for variations of 
the gradient V# m in the x direction (i.e., in the plane of 
the thin film) is large in comparison to the film thickness 
d, the resulting conservation of the momentum compo- 
nents k x and k y up to shifts of order 1/Zd <C l/d, together 
with energy conservation, constrains the possible values 
of the transverse momentum k z . If, in addition, the film 
thickness d is much smaller than the superconducting co- 
herence length £g = hvs/A, one then finds that the off- 
diagonal elements of the Andreev reflection amplitude 



r^(k y ;s;n' z s';n z s) are much smaller than the diagonal 
elements, so that one may consider r^(k y ;e;n' z s';n z s) 
to be diagonal in the mode indices n z and n' z . This is 
the case that we consider in the remainder of this section. 

As a first example, we consider the case that the mag- 
netization direction m has the spatial dependence (33) 
with 9 m a function of x only, so that the magnetization 
gradient is parallel to the HS interface. In that case one 
finds 



2v tx d 

x J dxr hc (x;k f e)e 2 * £X / hv ^, 

(52) 

where rh c (x; k|| ; e) is the Andreev reflection amplitude of 
Eq. (41), evaluated with the derivative d9 m /dx at posi- 
tion x. The prefactor v^ z /2v-\ x d in Eq. (52) expresses 
the geometric enhancement of the reflection amplitude 
from the coherent superposition of multiple reflections at 
the half-metal-superconductor interface. 32 The complex 
exponential factor accounts for the phase differences ac- 
quired by electrons and holes between these reflections. 
The in-plane wavevector k|| = sk x (0)e x + k y e y ; the small 
difference between the wavenumber k x (e) of the incom- 
ing electron and the wavenumber k x (— e) of the Andreev 
reflected hole is inconsequential for the calculation of r^ c 
in the Andreev approximation fcs£s S> 1- 

If the region in which the magnetization direction 
varies has a size smaller than the superconducting co- 
herence length £s , the x dependence of the complex ex- 
ponential e 2l£X / hv tx i n Eq. (52) can be neglected. Since 
rh e (a;;k||) is proportional to d9 m /dx, the integral over x 
gives the total change S9 m of the magnetization angle # m , 
and one finds 32 



r^(k y ;e;n' z s';n z s) 



\tf 



v^ z e 



i(4>-4>rL 



sk x (0)A 



+ 



2v tx dVA 2 - e 2 
v^ z Retti 



59„ 



4fc L ' {ktz + ^Vv^^u 



(53) 



In this limit, the Andreev reflection amplitude no longer 
depends on the size of the domain wall, nor on the pre- 
cise x dependence of the magnetization angle 9 m . In 
the opposite limit that the domain wall size Id is large 
in comparison to the superconducting coherence length, 
the reflection amplitude at e = is still given by Eq. (53) 
above, but Andreev reflection is suppressed for excitation 
energies e above the Thouless energy Ei d = Hv^ x /ld of 
the domain wall. The precise functional form of the sup- 
pression depends on the domain wall profile. A common 
domain wall profile is 46 



9m{x) = arctan[sinh(7ra;/Zd) 



(54) 



In this case the suppression at finite bias e — eV is by a 
factor 1/ cosh(eld/hvf x ). 
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Similar to Eq. (52) above, one finds the same result for 
variation of the magnetization along the direction normal 
to the superconducting interface, the z direction, as well 
as for a spin active interface. In these cases, rh c (x, y; k|| ) 
in Eq. (52) is the Andreev reflection amplitude of Eq. 
(40) and Eq. (44), respectively. 

VI. GREEN FUNCTIONS 

Most of the existing theoretical literature on the su- 
perconducting proximity effect in ferromagnets and half 
metals makes use of the Green function approach, with or 
without quasiclassical approximation. In the Green func- 
tion approach, the induced superconducting correlations 
in the half metal are characterized using the "anomalous 
Green function" . The symmetries of this Green function 
in the spin, orbital and, in particular, frequency domains 
arc used to classify the various forms of the proximity 
effect. 3 ' 6, 27,45 In this section, we investigate the relation 
between the scattering approach used here and the Green 
function approach. (The symmetry in the frequency do- 
main does not play an important role for the scatter- 
ing approach, because all information is encoded in "re- 
tarded" structures in the scattering approach.) 

The fundamental equation for the Green function 
Q(r, r';ifl) is the Gorkov equation 

(*fi-'W)0(r,r / ;tfi) = J(r-r / ), (55) 



where H is the Bogoliubov-de Gennes Hamiltonian of Eq. 
(9) and ft the Matsubara frequency. Like the Hamil- 
tonian H, the Green function Q has a 4 x 4 matrix 
structure, corresponding to the electron/hole and ma- 
jority/minority spin degrees of freedom. The anomalous 
Green function is the electron- hole (eh) block of Q. In a 
half metal, only the majority component is relevant, so 
one has 



F(r,r';ifl)=g etM (r,r';m). (56) 



We first calculate the anomalous Green function 
F(r, r'; ift) in the vicinity of the interface between a semi- 
infinite half metal and a superconductor. This is the 
geometry considered in Sec. IV. In order to calculate 
F(r, r';ift) we Fourier transform the Gorkov equation 
(55) to the coordinates x and y, and solve the remain- 
ing one-dimensional problem using the solutions of the 
Bogoliubov-de Gennes equation calculated in Sec. IV. 
For ft > and coordinates z, z' < (i.e., for coor- 
dinates inside the half metal), this procedure expresses 
the anomalous Green function F(r, r'; ift) in terms of the 
(outgoing) electron component of the exact hole-like re- 
tarded scattering state, and one finds 



i f e -ik z (z-z')+ik :L (x-x')+ik y (y-y')-\n\(\z\ + \z'\) 

F(r,r';ift) = dk x dk y r oh (k;ift), ft > 0, (57) 

(Z-rry J ihvf z 



where 



k z — ^Jk 2 — fc 2 — k 2 . (58) 



Similarly, for ft < 0, the anomalous Green function F(r, r'; ift) is found to be proportional to the (incoming) electron 
component of the exact hole-like advanced scattering state, 



F(r,r' ;J ft) = -jt-^t: I dk x dk y — r hc (k|| ; i|ft|)*, ft < 0. (59) 



I f e ik z (z-z')+ik x (x-x')+ik v (y-y')-\Q\(\z\ + \z'\) 

(2nfhj dkxdky 
These expressions can be cast in the form 



F(r,r';*ft) = — f dke ik '( r ~ r ')/(k, R; ifl)5(sk — Mht)> 
where R = \ (r + r'), Ek — h 2 k 2 /2mn, and 



(60) 



r oh (k||;z|ft|) if k z < and ft > 0, 
/(k,R;ift) = 2e- 2 l n ll R - e *l/ ft ^ x <j -r hc (kf if k z > and ft < 0, (61) 

otherwise. 

The function / may be identified with the anomalous Green function in the quasiclassical theory. 47 Note that / is 
nonzero only if the wavevector k points away from the superconductor interface if ft > (corresponding to a hole 
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moving towards the superconductor, retarded case), or if k points towards the superconductor if ft < (corresponding 
to an electron moving towards the superconductor, advanced case). 

Because of the delta function in Eq. (60), the function / is meaningful for wavevectors k with |k| = kf only. 
Following the literature, we analyze moments of /, taken with respect to its angular dependence on k, and determine 
whether these are even or odd functions of the Matsubara frequency ft. The moments are defined through the relation 

/(k,R;tfi) = ^ m (e k )/ im (R;iQ), (62) 

l.m 

where Yi m is a spherical harmonic and e k = k/fc the unit vector in the direction k. We restrict the discussion below 
to I = and 1 = 1. 

In order to determine the parity of the moments fi m , we note that 

r e h(k||,i|fi|) = rh e (-k||,i|fi|)*, 

see Sec. II. If the magnetization gradient is perpendicular to the half-metal-superconductor interface, rh (— kn , i\fl\) 
is an odd function of frequency, proportional to ft for small frequencies, but an even function of k||, see Eq. (40). 
Hence /oo is an odd function of ft, proportional to ft for small frequencies. 35,45 The moment /io is an even function 
of ft, proportional to |f2| for small ft. The functions fi m with m = ±1 are both zero. If the magnetization gradient is 
parallel to the half-metal-superconductor interface, r ne (— kii,i|fi|) is an even function of frequency, with a finite value 
for ft — > 0, but an odd function of k||, see Eq. (41). This implies that /oo and /io are both zero, whereas f\ m with 
m = ±1 are even functions of ft, with a finite nonzero value in the limit ft — »■ 0. 

The moments fi m calculated via the scattering approach obey the general symmetries imposed by the Pauli principle: 
For proximity effects induced by an s-wave, spin-singlet superconductor, the anomalous Green function fi m (fl) in a 
half metal is an odd function of ft for even I and an even function of it for odd I. 27 ' 45 Although foo(ft) — > for ft — > 
in the two cases discussed above, the requirement that the s-wave amplitude foo(fl) be an odd function of ft does 
not necessarily imply that /oo(^) always vanishes in this limit. Indeed, in junctions between a superconductor and a 
standard (not half-metallic) ferromagnct, it is known that the triplet component of foo(ft) approaches a finite value 
if ft — > 0. (To be precise, the triplet component of foo(ft) oc sign (ft) for small ft for a ferromagnet, in order to satisfy 
the antisymmetry constraint. 3 ' 26 ' 45 ) For junctions involving a half metal, an example of such a singular frequency 
dependence is given by the case of a lateral contact to a thin half-metallic film case, which we now discuss. 

As in Sec. V, we consider the case that the magnetization direction does not depend on y, and that the region 
of nonuniform magnetization is limited to the vicinity of x = 0. We assume that the inequalities described in the 
paragraph following Eq. (51) are obeyed, so that the thin-film reflection matrix r|f is diagonal in the transverse-mode 
indices n and n'. We calculate the anomalous Green function F(r, r';ifl) for the case that the coordinates r and r' 
are on the positive- x side of the region with non-uniform magnetization. Expanding the Green function in transverse 
modes, one then finds 

F(r,r';ifl) = ^-L, W dk„ sin ™ sin ^ e *ir('-0 /n ( k|| , R; m)d ( £k - m ), (63) 

n 

where 

( r c £(k y ;i\fi\;n+;n-) if k x > and ft > 0, 
/ n (k||,R;ifi) = 2e- 2 l°l( R - e ^/^ x I -rgj[(&„;i|fi|;n+;n-)* if k x < and ft < 0, (64) 

[ otherwise. 



The moments fi m with m + 1 odd vanish in the thin film 
geometry because they are odd in k z . Hence, we expand 
/ n (k|| , R; ift) in moments as 

/„(k||,R;in)= Y lm (e k )f lm (R;ifl), (65) 

l+m even 

where efe is the unit vector pointing in the direction 
k|| + nire z /d. In the thin-film geometry, the lowest mo- 
ment /oo is not only nonzero for magnetization gradients 



I 

perpendicular to the surface, but also for magnetization 
gradients parallel to the surface. In both cases /oo is an 
odd function of the frequency ft, but whereas /oo oc ft in 
the case of a magnetization gradient perpendicular to the 
interface, /oo is discontinuous, /oo oc sign(fi), in the case 
that the magnetization gradient is parallel to the surface. 

The fact that the s-wave amplitude /oo(w) may have 
a finite limit in the limit u> ~ > is a striking difference 
between the thin-film geometry and the regular geom- 
etry with a half-infinite half metal. Thus, our calcula- 
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tion identifies geometry as an important second factor 
in determining the symmetries of the anomalous Green's 
function — in addition to the Pauli principle. 

VII. CONCLUSION 

In this article we have calculated the Andreev reflection 
amplitudes of a half-metal-superconductor (HS) junction 
with a spatially nonuniform magnetization direction in 
the half metal. General symmetry considerations enforce 
that the Andreev reflection amplitude rh is zero at the 
Fermi level s = 0, except if inversion symmetry around 
the normal to the HS interface is broken. 29 ' 32 On the 
other hand, if that is the case, rh c is an odd function of 
the wavevector component k|| parallel to the interface. 
These general results were confirmed by explicit calcu- 
lations of rhc for magnetization gradients parallel and 
perpendicular to the interface. 

Impurity scattering has not been included in the calcu- 
lations presented here. As discussed in the introduction, 
this is not a serious shortcoming for the microscopic An- 
dreev reflection amplitude r c h (the amplitude for a single 
reflection of a quasiparticle incident on the HS interface) 
if the disorder is weak (mean free path I much larger 
than the Fermi wavelengths in the superconductor or in 
the half metal, and than the wavefunction decay length 
in the half metal), because Andreev reflection is a pro- 
cess that takes place on these microscopic length scales. 2 
Hence, the microscopic reflection amplitudes r c h of Sec. 
IV can be used as a valid starting point for a scattering 
theory of a disordered HS junction. The same situation 
occurs at normal-metal-superconductor (NS) interfaces, 
where the Andreev reflection amplitudes of a clean NS 
interface are combined with standard theoretical meth- 
ods for disordered normal metals in order to construct a 
theory of a disordered NS junction. 41 

However, disorder has a profound effect on the type 
and magnitude of the induced superconducting correla- 
tions in the half metal, which are mediated by the effect 
of multiple Andreev reflections. This includes the calcu- 
lation of the effective Andreev reflection amplitude r^f 
in the thin-film geometry, which represents the coherent 
effect of multiple Andreev reflections. In the absence of 
disorder, the wavevector component k|| parallel to the HS 
interface is conserved, and multiple Andreev reflections 
at the HS interface always add constructively. 32 With dis- 
order, k|| is no longer conserved. Because the sign of the 
Andreev reflection amplitudes r c h depends on the angle 
between k|| and the gradient of the magnetization angle, 
interference between multiple reflection events need no 
longer be constructive. In particular, the p-wave super- 
conducting correlations induced by a magnetization gra- 



dient parallel to the HS interface are strongly sensitive to 
disorder, and unable to survive for large distances away 
from the superconductor. 48 A detailed study of disorder 
effects in a lateral contact between a thin half-metallic 
film and a superconductor is left for a future publication. 

We close by noting that very recently half-metal- 
superconductor hybrid systems have garnered atten- 
tion as possible candidates for the creation of Ma- 
jorana fermion excitations, 49 which are considered 
promising candidates to implement topological quantum 
computing. 50,51 While the system considered here is too 
elementary to be useful for quantum computation, the 
fundamental ingredients present here — spin-flip scatter- 
ing, half-metallicity, and s-wave superconducting order 
— are precisely the same as those appearing in the pro- 
posals for Majorana fermion excitations. 49 Further anal- 
ysis of common and distinguishing features is thus desir- 
able. 
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Appendix A: Explicit calculation of the gauge 
transformation 

In this appendix we consider a general variation the 
magnetization direction, described by variations of both 
polar angles 9 m and 4> m . If both # m and </> m are position 
dependent, the transformation given in Eq. (35) leads to 
the gauge potential 

A = [(a 2 cos (j) m - (7i sin (f> m ) V6> m 

- ((Ti cos m sin # m + cr 2 sin </> m sin m ) V0 m 
-CT 3 (l-cos6» m )V0 m ]. (Al) 

In the calculation of Sec. IV only the terms propor- 
tional to V# m were included. There are three terms 
proportional to V0 m . Of these, the term proportional 
to 1T3 docs not contribute to spin-flip Andreev reflec- 
tion. Repeating the calculation of Sec. IV with the 
two remaining terms proportional to V^ m gives the re- 
sults (40) and (41) for gradients perpendicular and par- 
allel to the interface, respectively, with the replacement 
V# m — > V6* m + i sin6> m V(/> m . In the general case that 
the magnetization gradient has components parallel and 
perpendicular to the interface, the Andreev reflection am- 
plitudes of Eqs. (40) and (41) must be added. 
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